%函数：get_result(初始概率*100，概率开始增加的抽数，每次增加的概率*100）
%输出值：返回抽中的期望抽数和每一抽的期望概率
%此函数只适用于前N次概率不变，第K次开始概率均匀增加的抽卡模型

function [num,pro] = get_result(p_first,inflexion,p_add)
%初始概率，%概率转折点，%每次增加概率
p_add_num=0;%概率增加到1所需抽数
while p_first/100+(p_add_num+1)*p_add/100<1
    p_add_num=p_add_num+1;
end
A=ones(1,inflexion-1)*p_first/100;
B=linspace(p_add/100,p_add_num*p_add/100,p_add_num)+p_first/100;
P=[A,B,1];%第K次抽中的概率
Pi=ones(1,inflexion+p_add_num);a=1;q=0;
for k=1:(inflexion+p_add_num)
    Pi(k)=a*P(k);%恰好第K次抽中概率
    q=q+Pi(k)*k;%出货期望抽数
    a=a*(1-P(k));
end
plot(Pi);
format long
num = q;
pro = 1/q*100;
format default
